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ABSTRACT 

Adopting the Lindblom-Detweiler formalism for polar oscillations of neutron stars, we 
study the «;-mode oscillation and find that the Lagrangian change in pressure, mea- 
sured by the physical quantity X, is negligibly small. Based on this observation, we 
develop the generalised inverse-Cowling approximation (GICA) with the approxima- 
tion X = X' = 0, where X' is the derivative of X with respect to the circumferential 
radius, for w-mode oscillations of neutron stars. Under GICA, w-mode oscillations are 
described by a second-order differential system, which can yield accurate frequencies 
and damping rates of quasi-normal modes. 
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1 INTRODUCTION 

The search for gravitational waves emitted from stellar collapse has been the major goal of many relativists and astrophysicists 
since the early seventies of last century. At the beginning of a new millenium, several ground-based interferometric detectors 
(e.g. LIGO, VIR GO, GEO 600, TAMA 30 0) are currently endeavouring to detect gravitational waves from extraterrestrial 
objects (see, e.g., lHugheJ2003t iMasorfeooJ . and references therein). With the upgrade of these detectors and the availability of 
a space-based interferometric detector, LISA, which is expected to be launched at the beginning of the next decade, it is widely 
believed that gr avitational wave detection a nd hence astronomy in the gravitational wave window will become achievable in 
the near future (|Hughesll2003l . iMasonl liooil ) . In fact, there have been four science runs on LIGO and GEO 600 since 2002. 
Despite that so far no positive evidence for the existence of gravitational waves has been found, relevant obse rvational data 



were used to place limits on the ellipticity and the gravitational wave strain of a number of known pulsars (jAbbott et al 
20051 , 120071 ). 



One major kind of gravitational wave emitting astrophysical processes is certainly stellar gravitational collapse, often 
resulting in remnants of neutron stars or black holes, which ar e also possible grav itational wave emitters (for a review 
on gravitational waves emitted in gravitational collapse see, e.g., Fryer fc New 20031 . and references therein). For example, 
gravitational waves emitted in the formation of neutron stars in our galaxy via stellar core coll apses have been shown to 

Since such gravitational 



be detectable for sufficiently fierce asymmetric collapse ( Fryer et al. 20021 . 



Lindblom et al 



1998). 



waves are likely to carry the information about the internal structure of the neutron star involved, such as its mass, radius 
density distribution and in turn its equation of state (EOS), the characte ristics of gravitational wave signals emitted from 
undulating neutron stars have become a subje ct of much study (see, e.g.. 
19991 . iKokkotas et allboOll . iBenhar et aljliool ). 



Andersson fc Kokkotadll99d Il99gt benhar et al 



The spectrum of gravitational waves from compact stellar objects are commonly described in terms of quasi-normal modes 
(QNMs). Each of these QNMs is characteri z ed by a complex eigenfrequ ency uj — ui r + iuii and has a time dependence exp(iuit) 
! PresJl97ll . lLeaver|[l98rj . lChing et al.lll996l . lKokkotas fc Schmidtiri999l ). The QNM frequencies of gravitational waves emitted 
from a neutron star are generally model-sensitive and are expected to depend significantly on the EOS used in the stellar 
model. In spite of this fact, some universal behaviours in the frequency uj t and the damping time r = 1/ui of the /-mode and 
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the w-mode oscillations of non-rotating neutron stars can still be identified (jAndersson fc Kokkota Il998l . Bcnhar ct al. 1999), 
which have recently been explained by Tsui fc Leung (2005b). Based on such universal behaviours, it has been shown that 
the radius and th e mass of a neutron star can be inferred from the pulsation frequenci es of its fundamental fl uid / -mode and 
the fir st w-mode ( Andersson fc Kokkotasl 19981 . Benhar et al. 1999h . Most interestingly, Tsui fc Leung (2005a) and Tsui et al 



1 20061 ) showed that the EOS of a neutron star can be inferred from the eigenfrequencies of a few (say, e.g., 3) w-mode QNMs. 
These discoveries all underscore the physical significance of QNMs of neutron stars. On the other hand, it is important to note 
that detection of gravitation waves from neutron stars with s ufficient accuracy to infer th eir i nternal structure is a mission 
impossible in the short (or even medium) term. As estimated by Andersson fc Comer (2001) and Tsui et al. ( 20061 ). realisation 
of such schemes has to await the availability of more advanced gravitati onal-wave observatory lik e EUR O. 

Non-radial oscillations of relativistic neutron stars, first studied by Thorne fc Campolattarol 1 1967 ). can be analysed by 
decomposing the perturbed Einstein equations into spherical harmonics. Such oscillations are categorized into axial (toroidal) 
and polar (spheroidal) modes according to the parity of the harmonics. In axial-mode oscillations, which can be described 
by a single second-order wave equation IjChandrasekhar fc Ferraril 11991 ), fluid elements of a star a re merely spectators in 
the sense that they are not affected by gravitational waves generated IjThorne fc Campolattardll967l ). On the other hand, in 
polar-mode oscillations, fluid motion and gravitational wave are coupled together. Due to the complicated interplay between 
the matter and the metric, polar- mode oscillations are usuall y governed by a system of ordina ry differential equations (ODEs) 
involving several physical variables. In t he pioneering work of Thorne fc Campolattarol (|1967l ). polar oscillations are described 
by a fifth-order system of ODEs. Later, iLindblom fc Detweilerl (|1983l ) explicitly reduced the equation set into a fourth-order 
one, hereafter referred to as the LD formalism, where the oscillation of a star is described by four independent physical 
quantities Hi, K, W , and X. While Hi and K are measures of the metric, W and X are respectively proportional to the 
Lagrangian displacement in the radial dir ection and the pressure chang e of the fluid constituting the star (see Sect. 2 for more 
details). Besides the scheme proposed by Lindblom fc: Detweilei ll 19831), there are other schemes describing pol ar oscillations 
with two coupled second-order ODEs ( Ipser fc PriceHl99lT Price fc Ipser 199ll . Koiima 19921 . Allen et al. 19981 ). They are in 
fact equivalent to fourth-order systems. 

As in the Newtonian theory of stellar pulsations, polar oscillations can be classified into the fundamental (/) mode, the 
pressure (p) mode, and the gravity (g) mode. In addition to these, there is an extr a kind of mod e in the relativistic theory , 
namely the spacetime (w) mode, representing main ly perturbations in the metric ( Koiima 19881 . Kokkotas fc Schutz 1992 , 
Andersson et al. 19951 . Andersson fc Kokkotasl 19961 ) . 



Soon after the discovery of the p olar w-mode, it has been conjectured that the fl uid in a relativistic star will not be 

Aiming at understanding qualita- 



affected much by w-mode pulsations ( Kokkotas fc Schutzl 



1992, 



Andersson ct al 



1995). 



tively the physical nature of polar w-mode pulsations. I Andersson et alj l| 19961 ) suggested the Inverse-Cowling Approximation 
(ICA) for polar w-mode oscillations. In ICA the fluid motion is completely neglected from the outset (details of ICA will be 
discussed in Sect. 2 of the present paper) and the metric variables are governed by a second-order ODE system. The approx- 
imation can successfull y reproduce qualitative features of polar w-mode oscillations and clarify the nature of the oscillations 
I Andersson et al. 19961 ). However, as shown in Fig. 1, it fails to reproduce accurate values of the QNM frequencies. This is a 



sort of expected because the main concern of ICA was qualitative discussion on polar w-mode oscillations. 

The major objective of the present paper is to establish a quantitatively correct approximation for the polar w-mode, 
referred to as Generalised Inverse-Cowling Approximation (GICA) in the present paper, which is capable of locating polar w- 
mode with high accuracy (see Fig. 1). In our analysis we find that the fluid motion of a star is in fact non-negligible. However, 
X is so small that the eigenfrequency of polar w-mode can be accurately determined with the approximation X — X' — 
(hereafter a prime indicates differentiation with respect to circumferential radius r) by solving two coupled first-order ODEs. 
Thus, under GICA polar w-mode oscillations are describable with one single second-order ODE system. 

The organization of our paper is as follows. In Sect. 2, we briefly review the LD formalism. The ICA developed by 
Andersson et al. (1996) is discussed in Sect. 3 and we propose GICA in Sect. 4. We then conclude our paper in Sect. 5 with a 



short conclusion and discussion. Unless otherwise stated, geometrized units in which G = c = 1 are adopted in the following 
discussion. 



2 LD FORMALISM 

In spherical coordinates, (t,r,9,tfi), the geometry of spacetime around a non-rotating unperturbed neutron star is given by 
the line element: 

ds 2 = -e" (r) di 2 + e A(r) dr 2 + r 2 (d8 2 + sin 2 6»d^ 2 ). (1) 
The metric coefficient e A ' r ' is determined by the mass distribution function m(r), the mass inside circumferential radius r, 
e _ AW = x _2m(r) (2) 
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The other metric coefficient e"' r ' and t he mass distribution function m(r) can be ob tained from the solution to the Tolman- 
Oppenheimer-Volkoff (TOV) equations 1 Tolman 19391 . Oppenheimer fc Volkofl 19391 ): 



Av 2m + 8nr p 

dr r(r — 2m) ' 

dm 



dr 



(3) 
(4) 
(5) 



dp 1 Av 

d~r = ~2 (p + P) cV 
where p and p are the mass density and pressure, respectively. 

On the other hand, in the LD formalism for a pulsating relativistic star with radius R ( Lindblom fc Detweiler 19831 ). 
the perturbation in the metric, h^, can be expressed in terms of three physical quantities Ho, Hi and K and displayed in a 
matrix form as follows: 

/ -eV-Ho -iLup'Hi 





-iup l Hi -eV#o 
-p}r 2 K 





-fi l r 2 sin 2 8K J 



e ut Yi m (6, ifi), 



where to is eigenfrequency of the polar mode oscillation, I is the angular momentum quantum number, and 
[ 1 r^R. 

Correspondingly, the Lagrangian displacements of the fluid are expressed by: 

e = / u ! r- 1 e- A/ Vy !m (0, ¥ >)e iw * ) 

£ fl = —fi r~ 2 Vd$Yim(0,<p)e m ' i , 

C = -^{rsaxeyWd^Yi^e,^)^. 



(6) 



(7) 



(8) 
(9) 
(10) 



In other words, W and V together measure the spatial displacem ents of the fluid. 

Putting all these perturbations into the Einstein equations, Lindblom fc Detweiler ( 19831 ) showed that these physical 
variables are related by four first-order ODEs: 



Hi 



2me 2 \ . 

/ + H h 4nr e (p - p) 



K' = -H + 1-1(1 + l)Hi - 
r 2r 



2 + 1 



Hi + -e x [H + K - 16tt(p + p)V] 
87r(p + p)e A/2 1 



K - 



-W, 



W = 



l + i 



W + re 



A/2 



v/2 



IP 



-X 



1(1 + 1) 



V + -H + K 



X' = — X+(p + p)e 
r 



v/1 



~ - - ^ )H + i 



4n(p + p)e A//2 + o/e 



2 \/2-v 



2 T dr 



2 -„ , 1 1(1 + 1) 

tlo e H - 

2 r 



Hi + 



2 \2 rj 



(11) 
(12) 
(13) 

(14) 



where p(r), p(r) are respectively the density and the pressure at a radius r, 7 is the adiabatic index, defined by 
p + pdp 



7 = 



and 



P dp' 



(15) 



X = lo 2 ( p + p)e~ v/2 V - -p'e { "- x)/2 W + hp + p)e v/2 H . (16) 
r 2 

The physical significance of X is that it is a measure of Ap, the Lag r angia n ch ange in the fluid pressu r e. In fact, it can be 
readily shown from the theory developed bv lThorne fc Campolattard (1967) and lLindblom fc Detweilerl (1983) that: 

" //2 /i i XY !m (e,y.)e i " t . 



Ap = — e 

Besides, the metric perturbation Ho is expressed in terms of Hi, K and X 

Ho 

1 



1 * 
3m + -(I +2)(l - l)r + Airrp 



8nr 3 e-' /2 X 



— 1(1 + l)(m + 47rr 3 p) — u/Ve 



2 3 -(f + A) 



Hi 



— (I + 2)(l — l)r — uj 2 r 3 e " e A (m + 4-7rr 3 p)(3m — r + 4nr 3 f 

2 r 



K. 



(17) 



(18) 
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Eqs. - (|14|) . supplemented by the algebraic equations (|16[) and (1181) . t hen form a non- singular fourth-order system of 
ODEs with {Hi, K,W, X} being the four independent variables (jDetweiler fc Lindblomlll985T ). These four ODEs plus (i) the 
regularity condition at r = 0, (ii) the condition X — at r — R, and (iii) the outgoing bo undary condition at r = oo 
completely determine the eigenfrequency of polar mode pulsations ( Lindblom fc Detweilei 19831 ). 



3 INVERSE-COWLING APPROXIMATION 



Working in LP formalism and motivated by the conjecture that the fluid in a relat ivistic star is hardly excited by tu-mode 
pulsations ( Kokkotas fc Schutd 1992 . Andersson et al. 19951 ). Andersson et al. ( 19961 ) attempted to study the physical nature 
of polar w-mode oscillations by proposing the ICA with V = W = 0. Under such assumption, the motion of the fluid is 
completely ignored. To gauge the accuracy of the assumption V = W — 0, we show in Fig. 1 the I = 2 (always assumed for 
all numerical results in th is paper) polar w - mode QNM frequencies of (a) a polytropic star wi th p — up 2 and k = 100, which 
was used in the paper bv lAndersson et al. and (b) an APR1 star jAkmal et al.|[l998l ). which typifies other realistic 

neutron stars, and compare the numerical results obtained respectively from exact LD formalism, ICA with V = W — and 
GICA (to be developed in Sect. 4 of the present paper). It is obvious that in both cases the eigenfrequency obtained from 
the assumption V = W = deviates markedly from the exac t one and the agreement between the approximate and exact 
values is only qualitative. As noted by Andersson et al. ( 19961 ). the imaginary parts of QNM frequencies obtained from ICA 
are typically 20-30 percent smaller than the exact values and the relative error in the frequency spacing is about 10 percent. 
This is not really surprising because the original motivation of ICA was merely to verify the spacetime nature of polar w-mode 
oscillations. 

The inadequacy of the approximation V = W = can be clearly shown by comparing the contributions of various terms 
in (|16[) . whose physical meaning will be discussed at the end of this paper. In Fig. 2, the absolute values of ui 2 (p + p)e~"^ 2 V, 
p'e ( - I/ ~ x ^ 2 W/r, (p + p)e"/ 2 Ho/2 and X are compared for the a typical polar w-mo&e of the APR1 neutron star considered 
in Fig. lb. We note the term uj 2 (p + p)e~ l/ ^ 2 V , a measure of the fluid motion in the tangential direction, is comparable to 
the term (p + p)e u ^ 2 Hq/2, a perturbation in the metric. Besides, the term p' & u ~ x '' 2 W/r, measuring the motion in the radial 
direction, is not as small as expected. Judging from the results shown in Figs. 1 and 2, we deem that the approximation 
V — W — is unjustified. Instead, as suggested by the comparison shown in Fig. 2, the approximation X — X' — seems to 
be more appropriate. Similar behavi our also prevai l s in sta rs constructed with other realistic EOSs. We note that this point 
has been discussed and exploited by Tsui fc Leund (2005b) to explain the universality observed in the polar ™-mode. In the 
following section we will base on GICA with X = X' = to establish a second-order ODE system that is able to locate 
w-mode QNMs accurately. 



4 GENERALISED INVERSE-COWLING APPROXIMATION 

Under GICA where X = X' = 0, JUJl and (jT6j yield respectively 



= 

and 

V = 



-(---)#o + - 



2 _„ 1 1(1 + 1) 

no e H 

2 r 



a I i \ A/2 , 2 A/2-1/ J- 2 

47r(p + p)c ' + uj e ' r — 

2 dr 



Hi + --)K 
2 y 2 r' 



m + 4pnr 



-e e 
2a; 2 r 2 (—2m + r) lu 2 



W. 



(19) 



(20) 



Solving these two equations, we can express W in terms of the metric coefficients Hq, Hi and K: 



-5m + r - 12pnr 3 1 r . _ v 2 21 „ 

; ; K H — \ l(l + 1) + 2e r u \ H± 

2(r - 2m) 4 J 



c v l(l + l)(m + 47rpr 3 )] + r 2 Lu 2 (3m - r + 4irpr 3 ) 



H 



where 



D 



-4ne (p + p) + 



2(r - 2m)r 2 aj 2 

7m 2 — 47rr 4 (p — p + 47rp 2 r 2 ) — 4m (r + 27r pr 3 
r 2 (r — 2m) 2 



(21) 



(22) 
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Moreover, following directly from (|18[) and X = 0, Ho becomes 



H 



3m + ±(Z + 2) (I - l)r + 4p7rr 3 



-1(1 + 1) 



(m + 4npr 3 ) + c 



+ 



(Z - 1) (7 + 2) r - — (m + 4p7rr 3 ) (3m - r + 4p7rr 3 ) - r 3 uj 2 



K 



(23) 



As a result, V, W and Ho are all expressible in terms of Hi and K, which are governed by two first-order ODEs and 

To solve this second-order ODE system, we first discuss the boundary conditions at the center r = 0. Straightforward 
expansion about the center shows that Eqs. (|lip and (|12[) lead to: 

= -— Hi(0) + -[H o (0)+K(0)-16-K(p o + po)V(0)], 
r r 







1 



H o (0) 



1 



l(l + l)Hi(0) 



1 + 1 K(0) - 87r (P° + Po) 



r ' 2r ' r 
Comparing (I24p and (I25|) . we see that: 

W(0) = -IV(Q). 

In addition, (|23(l gives 

K(0) = H o (0), 

at the center. Putting (|26|) and (|27l) back into (|20|) . we have 



W(0), 



K(0) = ffo(0) 



8tt 



(3p + p) 



w(o). 



Finally, from (|24)| or (p5)) , 
1 



Hr(0) 



[2«i4:(0) + 167r(p + p)W / (0)]. 



(24) 
(25) 

(26) 
(27) 

(28) 
(29) 



Therefore, for a given W(0), which merely provides a scale of the solution, the values of Hi(0), K(0), Ho(0) and 1^(0) are 
all determined. Compared with the boundary conditions at the center, the boundary condition at the stellar surface r = R 
is much simpler. There are two independent solutions near the surface with arbitrarily assigned values of Hi(R) and K(R). 
Integrating and (|12[) outward from r — and inward from r — R, and matching the values of Hi, K at some intermediate 
value of r, say r = R/2, can numerically determine the solution of the second-order system up to a multiplicative constant. 

The QNM frequencies obtained from GICA are shown in Fig. 1 and Table 1, which agree nicely with the exact values. In 
particular, both the percentage errors in the frequency and the damping rate, E T and Ei, decrease with increasing frequency, 
reflecting the decrease in the Lagrangian change in pressure in high-frequency polar w-modes. Meanwhile, it is also worthy 
to note that even the ton mode, characterized by large damping rates, can be reproduced accurately. Hence, GICA proposed 
here, which assumes X = 0, clearly outperforms the ICA with the approximation V = W = 0. 



5 CONCLUSION AND DISCUSSION 

We have developed in the present paper an accurate scheme, GICA, to describe w-mode oscillations of neutron stars. Instead 
of completely ignoring the displacement of the fluid element of a star under consideration, we have shown that the term X, 
which is proportional to the Lagrangian change in the pressure, is negligible. Under GICA, w-mode oscillations of neutron 
stars are governed by two first-order ODE in the metric coefficients Hi and K. The displacements of matter, measured by V 
and W , are properly considered and found to be expressible in terms of Hi and K. Numerical results obtained directly from 
the approximation X — X' — clearly demonstrate the validity and accuracy of GICA (see Fig. 1 and Table 1). 

In addition to yielding reliable numerical values of QNM frequencies for w-mode oscillations, GICA proposed here also 
provides a proper physical picture for such oscillations. A star undergoing polar w-mode oscillations has non-negligible 
displacements in the tangential and radial directions. However, as shown in Fig. 2, such displacements are driven by the 
perturbed spacetime, instead of the excess pressure due to the motion of matter. 

It is well known that axial w-mode oscillations of compact stars, which are merely oscillat ions in the metric variables, 
are governed by a second-order Regge- Wheeler wave equation I Chandrasekhar fc Ferrari! 1991 ). Our finding that polar w- 



mode oscillations are adequately described by one single second-order ODE system thus puts polar and axial w-modes on 
an equal footing in the sense that both of them are governed by second-order ODE systems. As shown in the review by 
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Kokkotas fc Schmidtl (|l999h . QNMs o f these two kinds of w-m ode lie asymptot i cally on a common smooth curve in the 
complex frequency plane. In addition, Koiima et al. (1995) and Andersson et al. ( 19961 ) also noted that for incompressible 
ultra compact stars the oscillation frequencies of these two iu-modes agree to within one percent. GICA proposed here may 
be the first step in the long process of seeking quantitative interpretation of such interesting asymptotic behaviour and the 
similarity between two kinds of w-mode. 

We have adopted LD formalism to describe polar oscill ations in the present paper. In fact, there are other alternative 



We have adopted Li) lormahsm to describe polar osculations in the present paper. In tact, there are othe r alternative 
form alisms for polar oscillations of relativistic sta rs (see, e.g., Ipser fc Price! 1991 . Koiima 1992 . Allen et al. 199§l . Nagar et al 



2004). In particular, in the formalism proposed by Allen et al. 1 1998h . and later used by Nagar et al.l ( 2004 ) to study accretion- 



driven gravitational waves from compact stars, such oscillations are described by two coupled second-order wave equations 
in the metric variables plus a Hamiltonian constraint yielding the Eulerian change in the fluid pressure. The approach is 
also equivalent to a fourth-order differential equation system with the tortoise ordinate being the independent variable. The 
formalism is evidently more amenable to decoupling metric and matter variables. It is worthy to develop quantitatively 
accurate approximation scheme for polar w-mode oscillations based on such formalism. Currently we are working along this 
direction. 

Lastly, it is worthy to note the physical meaning of 1)160 . which can be rewritten as follows: 



AtV(p + p) VY lm (9, tp)&* = -e v Sp - 
where 



p l (p + p)e l, H 



Yi m (0,(p)e 



8p = Ap - 



Ap- 



V l p'e 



-V2iy 



Yi m (6, ip)e lu 



(30) 



(31) 



is the Eulerian change in p ressure. Each term in (|30)l has an obvious analogue in Newtonian theory of stellar pulsation (see, e.g. 
Eq. (17.21b) in[Cox|[l980j). The term fj, l u> 2 (p + p)VYim(6, ¥')e 1 " t in the LHS is a measure of inertia x acceleration, whereas the 



two terms in the RHS, —e v 8p and — fj 1 (p + p)e" HoYi m (8 , ip)e lult /2, correspond to the pressure and gravity forces, respectively. 
Thus, (|30)1 can be understood as the extension of application of Newton's second law to general relativistic theory of stellar 
pulsation. Following directly from Fig. 2 and 1)30)1 , the net driving forces for io-mode pulsation is the sum of gravity and the 
Eulerian change in pressure. 

To compare w-mode with / and p modes, we show in Fig. 3 and Fig. 4 the corresponding physical quantities appearing 
in 1)16)1 for the /-mode and the first p-mode of the same star. It is interesting to note that for /-mode oscillations the major 
driving force is the Eulerian change in pressure 8p and the gravity in fact partially cancels the effect of Sp. On the other hand, 
the major driving force for p-mode oscillations is evidently due to the Lagrangian change in pressure Ap. We expect that 
these findings can cast light on eigen-mode analysis for pulsations of relativistic stars and assist in developing appropriate 
approximation schemes under different circumstances. Relevant study is now underway and will be reported elsewhere in due 
course. 
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Table 1. Polar «j-mode QNM frequencies of the two typical stars considered in Fig. 1 are obtained from GICA and exact LD formalism, 
denoted respectively by w and w, and tabulated here. The percentage errors in the frequency and damping rate, E T = \w T /ui T — 1| X 100% 
and E[ = \u>i/u>i — 1| X 100%, are also included. 



Polytropic star 





Muji 






E T 


Ei 


0.0360 


0.7045 


0.0385 


0.7064 


6.5 


0.27 


0.3377 


0.3804 


0.3427 


0.3781 


1.5 


0.6 


0.5027 


0.2593 


0.5029 


0.2645 


0.022 


2.0 


0.8741 


0.3637 


0.8752 


0.3694 


0.12 


1.5 


1.2235 


0.4193 


1.2241 


0.4250 


0.05 


1.3 


1.5683 


0.4596 


1.5687 


0.4652 


0.026 


1.2 


1.9111 


0.4915 


1.9114 


0.4971 


0.016 


1.1 


2.2528 


0.5180 


2.2531 


0.5235 


0.011 


1.1 


2.5938 


0.5407 


2.5940 


0.5462 


0.0077 


1.0 


2.9343 


0.5606 


2.9344 


0.5660 


0.0059 


0.97 


3.2743 


0.5783 


3.2745 


0.5838 


0.0051 


0.94 


3.6141 


0.5943 


3.6143 


0.5998 


0.0046 


0.91 


3.9536 


0.6090 


3.9538 


0.6144 


0.0041 


0.88 


APR1 EOS 


MCj t 


Mu>i 


Muj r 


Mu>i 


E T 


Ei 


0.2963 


0.3551 


0.2983 


0.3533 


0.69 


0.51 


0.5057 


0.2924 


0.5073 


0.2963 


0.3 


1.3 


0.9247 


0.3777 


0.9266 


0.3818 


0.21 


1.1 


1.3138 


0.4307 


1.3155 


0.4353 


0.13 


1.1 


1.6958 


0.4700 


1.6973 


0.4749 


0.086 


1.0 


2.0754 


0.5001 


2.0767 


0.5051 


0.063 


1.0 


2.4551 


0.5238 


2.4564 


0.5287 


0.054 


0.93 


2.8361 


0.5431 


2.8376 


0.5478 


0.053 


0.85 


3.2188 


0.5594 


3.2205 


0.5637 


0.055 


0.77 


3.6028 


0.5733 


3.6049 


0.5771 


0.058 


0.66 
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Figure 1. Polar ui-mode QNMs of (a) a polytropic star with 7 = 2 and compactness = 0.211; and (b) an APR1 star with compactness = 
0.2 are obtained from exact LD formalism, ICA proposed previously and GICA developed here, respectively. Here M is the mass of the 
star. 




Figure 2. The absolute values of various terms in {TfiJ, u) 2 (p + p)c V / 2 V, p'ef'" ^'' 2 W/r, (p + p)e"/ 2 Ho/2 and X, are compared for 
the least damped polar ui-mode (I = 2) of the APR! neutron star considered in Fig. lb. 
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Figure 3. The absolute values of various terms in H16I I, u) 2 (p + p)e "/ 2 V, p'e(" x ^ 2 W/r, (p + p)e v ^ 2 and X. are compared for 
the /-mode (I = 2) of the APR1 neutron star considered in Fig. lb. 




Figure 4. The absolute values of various terms in QU, cu 2 (p + p)e _I/ / 2 V, p' e^'^^W /r , (p + p)e"/ 2 H /2 and X, are compared for 
the first p-mode (/ = 2) of the APR1 neutron star considered in Fig. lb. 



